Whole-body metabolic modelling reveals microbiome and genomic interactions on reduced urine formate levels in Alzheimer's disease.

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimer's disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilised whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalised whole-body metabolic models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

Replication Protein A, the Main Eukaryotic Single-Stranded DNA Binding Protein, a Focal Point in Cellular DNA Metabolism.

Replication protein A (RPA) is a heterotrimeric protein complex and the main single-stranded DNA (ssDNA)-binding protein in eukaryotes. RPA has key functions in most of the DNA-associated metabolic pathways and DNA damage signalling. Its high affinity for ssDNA helps to stabilise ssDNA structures and protect the DNA sequence from nuclease attacks. RPA consists of multiple DNA-binding domains which are oligonucleotide/oligosaccharide-binding (OB)-folds that are responsible for DNA binding and interactions with proteins. These RPA-ssDNA and RPA-protein interactions are crucial for DNA replication, DNA repair, DNA damage signalling, and the conservation of the genetic information of cells. Proteins such as ATR use RPA to locate to regions of DNA damage for DNA damage signalling. The recruitment of nucleases and DNA exchange factors to sites of double-strand breaks are also an important RPA function to ensure effective DNA recombination to correct these DNA lesions. Due to its high affinity to ssDNA, RPA's removal from ssDNA is of central importance to allow these metabolic pathways to proceed, and processes to exchange RPA against downstream factors are established in all eukaryotes. These faceted and multi-layered functions of RPA as well as its role in a variety of human diseases will be discussed.

A toolbox of machine learning software to support microbiome analysis.

The human microbiome has become an area of intense research due to its potential impact on human health. However, the analysis and interpretation of this data have proven to be challenging due to its complexity and high dimensionality. Machine learning (ML) algorithms can process vast amounts of data to uncover informative patterns and relationships within the data, even with limited prior knowledge. Therefore, there has been a rapid growth in the development of software specifically designed for the analysis and interpretation of microbiome data using ML techniques. These software incorporate a wide range of ML algorithms for clustering, classification, regression, or feature selection, to identify microbial patterns and relationships within the data and generate predictive models. This rapid development with a constant need for new developments and integration of new features require efforts into compile, catalog and classify these tools to create infrastructures and services with easy, transparent, and trustable standards. Here we review the state-of-the-art for ML tools applied in human microbiome studies, performed as part of the COST Action ML4Microbiome activities. This scoping review focuses on ML based software and framework resources currently available for the analysis of microbiome data in humans. The aim is to support microbiologists and biomedical scientists to go deeper into specialized resources that integrate ML techniques and facilitate future benchmarking to create standards for the analysis of microbiome data. The software resources are organized based on the type of analysis they were developed for and the ML techniques they implement. A description of each software with examples of usage is provided including comments about pitfalls and lacks in the usage of software based on ML methods in relation to microbiome data that need to be considered by developers and users. This review represents an extensive compilation to date, offering valuable insights and guidance for researchers interested in leveraging ML approaches for microbiome analysis.

Causal inference on microbiome-metabolome relations in observational host-microbiome data via in silico in vivo association pattern analyses.

Understanding the effects of the microbiome on the host's metabolism is core to enlightening the role of the microbiome in health and disease. Herein, we develop the paradigm of in silico in vivo association pattern analyses, combining microbiome metabolome association studies with in silico constraint-based community modeling. Via theoretical dissection of confounding and causal paths, we show that in silico in vivo association pattern analyses allow for causal inference on microbiome-metabolome relations in observational data. We justify the corresponding theoretical criterion by structural equation modeling of host-microbiome systems, integrating deterministic microbiome community modeling into population statistics approaches. We show the feasibility of our approach on a published multi-omics dataset (n = 347), demonstrating causal microbiome-metabolite relations for 26 out of 54 fecal metabolites. In summary, we generate a promising approach for causal inference in metabolic host-microbiome interactions by integrating hypothesis-free screening association studies with knowledge-based in silico modeling.

APOLLO: A genome-scale metabolic reconstruction resource of 247,092 diverse human microbes spanning multiple continents, age groups, and body sites.

Computational modelling of microbiome metabolism has proved instrumental to catalyse our understanding of diet-host-microbiome-disease interactions through the interrogation of mechanistic, strain- and molecule-resolved metabolic models. We present APOLLO, a resource of 247,092 human microbial genome-scale metabolic reconstructions spanning 19 phyla and accounting for microbial genomes from 34 countries, all age groups, and five body sites. We explored the metabolic potential of the reconstructed strains and developed a machine learning classifier able to predict with high accuracy the taxonomic strain assignments. We also built 14,451 sample-specific microbial community models, which could be stratified by body site, age, and disease states. Finally, we predicted faecal metabolites enriched or depleted in gut microbiomes of people with Crohn's disease, Parkinson disease, and undernourished children. APOLLO is compatible with the human whole-body models, and thus, provide unprecedented opportunities for systems-level modelling of personalised host-microbiome co-metabolism. APOLLO will be freely available under https://www.vmh.life/.

Whole-body modelling reveals microbiome and genomic interactions on reduced urine formate levels in Alzheimer's disease.

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimer's disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilized whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalized models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

Cardinality optimization in constraint-based modelling: application to human metabolism.

MOTIVATION: Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution. RESULTS: By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction. AVAILABILITY AND IMPLEMENTATION: Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.

The Effects of Hospitalisation on the Serum Metabolome in COVID-19 Patients.

COVID-19, a systemic multi-organ disease resulting from infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is known to result in a wide array of disease outcomes, ranging from asymptomatic to fatal. Despite persistent progress, there is a continued need for more accurate determinants of disease outcomes, including post-acute symptoms after COVID-19. In this study, we characterised the serum metabolomic changes due to hospitalisation and COVID-19 disease progression by mapping the serum metabolomic trajectories of 71 newly hospitalised moderate and severe patients in their first week after hospitalisation. These 71 patients were spread out over three hospitals in Switzerland, enabling us to meta-analyse the metabolomic trajectories and filter consistently changing metabolites. Additionally, we investigated differential metabolite-metabolite trajectories between fatal, severe, and moderate disease outcomes to find prognostic markers of disease severity. We found drastic changes in serum metabolite concentrations for 448 out of the 901 metabolites. These results included markers of hospitalisation, such as environmental exposures, dietary changes, and altered drug administration, but also possible markers of physiological functioning, including carboxyethyl-GABA and fibrinopeptides, which might be prognostic for worsening lung injury. Possible markers of disease progression included altered urea cycle metabolites and metabolites of the tricarboxylic acid (TCA) cycle, indicating a SARS-CoV-2-induced reprogramming of the host metabolism. Glycerophosphorylcholine was identified as a potential marker of disease severity. Taken together, this study describes the metabolome-wide changes due to hospitalisation and COVID-19 disease progression. Moreover, we propose a wide range of novel potential biomarkers for monitoring COVID-19 disease course, both dependent and independent of the severity.

Circulating metabolites modulated by diet are associated with depression.

Metabolome reflects the interplay of genome and exposome at molecular level and thus can provide deep insights into the pathogenesis of a complex disease like major depression. To identify metabolites associated with depression we performed a metabolome-wide association analysis in 13,596 participants from five European population-based cohorts characterized for depression, and circulating metabolites using ultra high-performance liquid chromatography/tandem accurate mass spectrometry (UHPLC/MS/MS) based Metabolon platform. We tested 806 metabolites covering a wide range of biochemical processes including those involved in lipid, amino-acid, energy, carbohydrate, xenobiotic and vitamin metabolism for their association with depression. In a conservative model adjusting for life style factors and cardiovascular and antidepressant medication use we identified 8 metabolites, including 6 novel, significantly associated with depression. In individuals with depression, increased levels of retinol (vitamin A), 1-palmitoyl-2-palmitoleoyl-GPC (16:0/16:1) (lecithin) and mannitol/sorbitol and lower levels of hippurate, 4-hydroxycoumarin, 2-aminooctanoate (alpha-aminocaprylic acid), 10-undecenoate (11:1n1) (undecylenic acid), 1-linoleoyl-GPA (18:2) (lysophosphatidic acid; LPA 18:2) are observed. These metabolites are either directly food derived or are products of host and gut microbial metabolism of food-derived products. Our Mendelian randomization analysis suggests that low hippurate levels may be in the causal pathway leading towards depression. Our findings highlight putative actionable targets for depression prevention that are easily modifiable through diet interventions.

Longitudinal flux balance analyses of a patient with episodic colonic inflammation reveals microbiome metabolic dynamics.

We report the first use of constraint-based microbial community modeling on a single individual with episodic inflammation of the gastrointestinal tract, who has a well documented set of colonic inflammatory biomarkers, as well as metagenomically-sequenced fecal time series covering seven dates over 16 months. Between the first two time steps the individual was treated with both steroids and antibiotics. Our methodology enabled us to identify numerous time-correlated microbial species and metabolites. We found that the individual's dynamical microbial ecology in the disease state led to time-varying in silico overproduction, compared to healthy controls, of more than 24 biologically important metabolites, including methane, thiamine, formaldehyde, trimethylamine N-oxide, folic acid, serotonin, histamine, and tryptamine. The microbe-metabolite contribution analysis revealed that some Dialister species changed metabolic pathways according to the inflammation phases. At the first time point, characterized by the highest levels of serum (complex reactive protein) and fecal (calprotectin) inflammation biomarkers, they produced L-serine or formate. The production of the compounds, through a cascade effect, was mediated by the interaction with pathogenic Escherichia coli strains and Desulfovibrio piger. We integrated the microbial community metabolic models of each time point with a male whole-body, organ-resolved model of human metabolism to track the metabolic consequences of dysbiosis at different body sites. The presence of D. piger in the gut microbiome influenced the sulfur metabolism with a domino effect affecting the liver. These results revealed large longitudinal variations in an individual's gut microbiome ecology and metabolite production, potentially impacting other organs in the body. Future simulations with more time points from an individual could permit us to assess how external drivers, such as diet change or medical interventions, drive microbial community dynamics.

Systems biology elucidates the distinctive metabolic niche filled by the human gut microbe Eggerthella lenta.

Human gut bacteria perform diverse metabolic functions with consequences for host health. The prevalent and disease-linked Actinobacterium Eggerthella lenta performs several unusual chemical transformations, but it does not metabolize sugars and its core growth strategy remains unclear. To obtain a comprehensive view of the metabolic network of E. lenta, we generated several complementary resources: defined culture media, metabolomics profiles of strain isolates, and a curated genome-scale metabolic reconstruction. Stable isotope-resolved metabolomics revealed that E. lenta uses acetate as a key carbon source while catabolizing arginine to generate ATP, traits which could be recapitulated in silico by our updated metabolic model. We compared these in vitro findings with metabolite shifts observed in E. lenta-colonized gnotobiotic mice, identifying shared signatures across environments and highlighting catabolism of the host signaling metabolite agmatine as an alternative energy pathway. Together, our results elucidate a distinctive metabolic niche filled by E. lenta in the gut ecosystem. Our culture media formulations, atlas of metabolomics data, and genome-scale metabolic reconstructions form a freely available collection of resources to support further study of the biology of this prevalent gut bacterium.

A nutrition algorithm to optimize feed and medium composition using genome-scale metabolic models.

The optimization of animal feeds and cell culture media are problems of interest to a wide range of industries and scientific disciplines. Both problems are dictated by the properties of an organism's metabolism. However, due to the tremendous complexity of metabolic systems, it can be difficult to predict how metabolism will respond to changes in nutrient availability. A common tool used to capture the complexity of metabolism in a computational framework is a genome-scale metabolic model (GEM). GEMs are useful for predicting the fluxes of reactions within an organism's metabolism. To optimize feed or media, in silico experiments can be performed with GEMs by systematically varying nutritional constraints and predicting metabolic activity. In this way, the influence of various nutritional changes on metabolic outcomes can be evaluated. However, this methodology does not guarantee an optimal solution. Here, we develop a nutrition algorithm that utilizes linear programming to search the entire flux solution space of possible dietary intervention strategies to identify the most efficient changes to nutrition for a desirable metabolic outcome. We illustrate the utility of the nutrition algorithm on GEMs of Atlantic salmon (Salmo salar) and Chinese hamster ovary (CHO) cell metabolism and find that the nutrition algorithm makes predictions that not only align with experimental findings but reveal new insights into promising feeding strategies. We show that the nutrition algorithm is highly versatile and customizable to meet the user's needs. For instance, we demonstrate that the nutrition algorithm can be used to predict feed/media compositions that maximize profit margins. While the nutrition algorithm can be used to define an optimal feed/medium ab initio, it can also identify minimal changes to be made to an existing feed/medium to drive the largest metabolic shift. Moreover, the nutrition algorithm can target multiple metabolic pathways simultaneously with only a marginal increase in computational expense. While the nutrition algorithm has its limitations, we believe that this tool can be leveraged in a broad range of biotechnological applications to enhance the feed/medium optimization process.

Preterm birth is associated with xenobiotics and predicted by the vaginal metabolome.

Spontaneous preterm birth (sPTB) is a leading cause of maternal and neonatal morbidity and mortality, yet its prevention and early risk stratification are limited. Previous investigations have suggested that vaginal microbes and metabolites may be implicated in sPTB. Here we performed untargeted metabolomics on 232 second-trimester vaginal samples, 80 from pregnancies ending preterm. We find multiple associations between vaginal metabolites and subsequent preterm birth, and propose that several of these metabolites, including diethanolamine and ethyl glucoside, are exogenous. We observe associations between the metabolome and microbiome profiles previously obtained using 16S ribosomal RNA amplicon sequencing, including correlations between bacteria considered suboptimal, such as Gardnerella vaginalis, and metabolites enriched in term pregnancies, such as tyramine. We investigate these associations using metabolic models. We use machine learning models to predict sPTB risk from metabolite levels, weeks to months before birth, with good accuracy (area under receiver operating characteristic curve of 0.78). These models, which we validate using two external cohorts, are more accurate than microbiome-based and maternal covariates-based models (area under receiver operating characteristic curve of 0.55-0.59). Our results demonstrate the potential of vaginal metabolites as early biomarkers of sPTB and highlight exogenous exposures as potential risk factors for prematurity.

Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine.

The human microbiome influences the efficacy and safety of a wide variety of commonly prescribed drugs. Designing precision medicine approaches that incorporate microbial metabolism would require strain- and molecule-resolved, scalable computational modeling. Here, we extend our previous resource of genome-scale metabolic reconstructions of human gut microorganisms with a greatly expanded version. AGORA2 (assembly of gut organisms through reconstruction and analysis, version 2) accounts for 7,302 strains, includes strain-resolved drug degradation and biotransformation capabilities for 98 drugs, and was extensively curated based on comparative genomics and literature searches. The microbial reconstructions performed very well against three independently assembled experimental datasets with an accuracy of 0.72 to 0.84, surpassing other reconstruction resources and predicted known microbial drug transformations with an accuracy of 0.81. We demonstrate that AGORA2 enables personalized, strain-resolved modeling by predicting the drug conversion potential of the gut microbiomes from 616 patients with colorectal cancer and controls, which greatly varied between individuals and correlated with age, sex, body mass index and disease stages. AGORA2 serves as a knowledge base for the human microbiome and paves the way to personalized, predictive analysis of host-microbiome metabolic interactions.

MetaboAnnotator: an efficient toolbox to annotate metabolites in genome-scale metabolic reconstructions.

MOTIVATION: Genome-scale metabolic reconstructions have been assembled for thousands of organisms using a wide range of tools. However, metabolite annotations, required to compare and link metabolites between reconstructions, remain incomplete. Here, we aim to further extend metabolite annotation coverage using various databases and chemoinformatic approaches. RESULTS: We developed a COBRA toolbox extension, deemed MetaboAnnotator, which facilitates the comprehensive annotation of metabolites with database independent and dependent identifiers, obtains molecular structure files, and calculates metabolite formula and charge at pH 7.2. The resulting metabolite annotations allow for subsequent cross-mapping between reconstructions and mapping of, e.g., metabolomic data. AVAILABILITY AND IMPLEMENTATION: MetaboAnnotator and tutorials are freely available at https://github.com/opencobra. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Whole-body metabolic modelling predicts isoleucine dependency of SARS-CoV-2 replication.

We aimed at investigating host-virus co-metabolism during SARS-CoV-2 infection. Therefore, we extended comprehensive sex-specific, whole-body organ resolved models of human metabolism with the necessary reactions to replicate SARS-CoV-2 in the lung as well as selected peripheral organs. Using this comprehensive host-virus model, we obtained the following key results: 1. The predicted maximal possible virus shedding rate was limited by isoleucine availability. 2. The supported initial viral load depended on the increase in CD4+ T-cells, consistent with the literature. 3. During viral infection, the whole-body metabolism changed including the blood metabolome, which agreed well with metabolomic studies from COVID-19 patients and healthy controls. 4. The virus shedding rate could be reduced by either inhibition of the guanylate kinase 1 or availability of amino acids, e.g., in the diet. 5. The virus variants differed in their maximal possible virus shedding rates, which could be inversely linked to isoleucine occurrences in the sequences. Taken together, this study presents the metabolic crosstalk between host and virus and emphasises the role of amino acid metabolism during SARS-CoV-2 infection, in particular of isoleucine. As such, it provides an example of how computational modelling can complement more canonical approaches to gain insight into host-virus crosstalk and to identify potential therapeutic strategies.

NMR Metabolomics Reveal Urine Markers of Microbiome Diversity and Identify Benzoate Metabolism as a Mediator between High Microbial Alpha Diversity and Metabolic Health.

Microbial metabolites measured using NMR may serve as markers for physiological or pathological host-microbe interactions and possibly mediate the beneficial effects of microbiome diversity. Yet, comprehensive analyses of gut microbiome data and the urine NMR metabolome from large general population cohorts are missing. Here, we report the associations between gut microbiota abundances or metrics of alpha diversity, quantified from stool samples using 16S rRNA gene sequencing, with targeted urine NMR metabolites measures from 951 participants of the Study of Health in Pomerania (SHIP). We detected significant genus-metabolite associations for hippurate, succinate, indoxyl sulfate, and formate. Moreover, while replicating the previously reported association between hippurate and measures of alpha diversity, we identified formate and 4-hydroxyphenylacetate as novel markers of gut microbiome alpha diversity. Next, we predicted the urinary concentrations of each metabolite using genus abundances via an elastic net regression methodology. We found profound associations of the microbiome-based hippurate prediction score with markers of liver injury, inflammation, and metabolic health. Moreover, the microbiome-based prediction score for hippurate completely mediated the clinical association pattern of microbial diversity, hinting at a role of benzoate metabolism underlying the positive associations between high alpha diversity and healthy states. In conclusion, large-scale NMR urine metabolomics delivered novel insights into metabolic host-microbiome interactions, identifying pathways of benzoate metabolism as relevant candidates mediating the beneficial health effects of high microbial alpha diversity.

The gut microbial metabolite formate exacerbates colorectal cancer progression.

The gut microbiome is a key player in the immunomodulatory and protumorigenic microenvironment during colorectal cancer (CRC), as different gut-derived bacteria can induce tumour growth. However, the crosstalk between the gut microbiome and the host in relation to tumour cell metabolism remains largely unexplored. Here we show that formate, a metabolite produced by the CRC-associated bacterium Fusobacterium nucleatum, promotes CRC development. We describe molecular signatures linking CRC phenotypes with Fusobacterium abundance. Cocultures of F. nucleatum with patient-derived CRC cells display protumorigenic effects, along with a metabolic shift towards increased formate secretion and cancer glutamine metabolism. We further show that microbiome-derived formate drives CRC tumour invasion by triggering AhR signalling, while increasing cancer stemness. Finally, F. nucleatum or formate treatment in mice leads to increased tumour incidence or size, and Th17 cell expansion, which can favour proinflammatory profiles. Moving beyond observational studies, we identify formate as a gut-derived oncometabolite that is relevant for CRC progression.

Microbiome Modelling Toolbox 2.0: efficient, tractable modelling of microbiome communities.

MOTIVATION: Constraint-Based Reconstruction and Analysis (COBRA) is a widely used approach for the interrogation and stratification of microbiome samples, yet applications to large-scale cohorts are hampered by limited scalability and efficiency of simulations. RESULTS: We substantially improved the computation speed and scalability of a previous implementation for the construction and interrogation of personalized constraint-based microbiome models as well as implemented additional functionalities for analysis and visualization. AVAILABILITY AND IMPLEMENTATION: Microbiome Modelling Toolbox and tutorials are freely available as part of the COBRA Toolbox at https://git.io/microbiomeModelingToolbox.